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Extensive computations of ground state energies of the Edwards- Anderson spin glass on bond- 
diluted, hypercubic lattices are conducted in dimensions d — 3, ... ,7. Results are presented for 
bond-densities exactly at the percolation threshold, p — p c , and deep within the glassy regime, 
p > p c , where finding ground-states is one of the hardest combinatorial optimization problems. 
Finite-size corrections of the form \/N" are shown to be consistent throughout with the prediction 
id = 1 — y/d, where y refers to the "stiffness" exponent that controls the formation of domain wall 
excitations at low temperatures. At p — p c , an extrapolation for d — > oo appears to match our 
mean-field results for these corrections. In the glassy phase, however, ui does not approach its 
anticipated mean-field value of 2/3, obtained from simulations of the Sherrington-Kirkpatrick spin 
glass on an TV-clique graph. Instead, the value of oj reached at the upper critical dimension matches 
another type of mean-field spin glass models, namely those on sparse random networks of regular 
degree called Bethe lattices. 
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I. INTRODUCTION 

The relevance of mean-field predictions based on 
the Sherrington-Kirkpatrick model (SK)i for the finite- 
dimensional Ising spin-glass introduced by Edwards 
and Anderson (EA) 2 - has been an issue of extensive 
discussions Qualitative predictions of mean-field cal- 
culations often are taken for granted in non-disordered 
systems. Yet, many questions have been raised££&ii~— 
about those predictions for the phase diagram of EA ob- 
tained from SK, as solved with replica symmetry break- 
ing (RSB) by Parisi^ii Whatever the true nature of the 
broken symmetry in a real Ising glass may turn out to be, 
one would expect tantalizing insights into the properties 
of energy landscapes of disordered systems generally^ by 
understanding how any mean-field RSB-solution morphs 
into its finite-dimensional counterpart across the upper 
critical dimension, here, d u = 6. Unfortunately, a direct 
study of this connection with an expansion in e = d u — d 
is beset with technical difficulties^ and results remain 
hard to interpret. Numerical approaches have renewed 
interest in this connection, based on simulations of Id 
long-range model a 12 ! 16 " — or of bond-diluted hypercubic 
lattices^ - — which we adapt for our study here. Prior 
to the development of this method, scaling studies on 
ground states beyond d — 4 (L < 7)22- were unheard of. 

In this Letter we investigate finite-size corrections 
(FSC) in EA on dilute lattices at some bond-density 
p and sizes N = L d in up to d = 7 dimensions and 
probe their connection with similar studies in mean-field 
models, such as SK or sparse regular random graphs fre- 
quently called Bethe lattices (BL)<2ir— For lattices at 
the bond-percolation threshold, the exponent u> adheres 
to the predicted scaling relation^ Eq. © below, and 
connects smoothly with its mean-field counterpart above 
di . Most of our computational effort is directed at show- 
ing that, in the glassy phase, id is equally consistent with 
the same scaling relation. Although widely conjectured, a 
rigorous argument for that relation at T = when T c > 



is lacking. The extrapolation to its apparent mean-field 
limit is surprising: it does not match up (nor cross over—) 
with corrections found numerically for either, SK— ~— or 
BL with ±J-bondSf2& Only FSC for BL with Gaussian 
bonds^ 7 - seem to provide a plausible mean-field limit for 
u> in EA. This evidence for a disconnect in the FSC be- 
tween EA and SK beyond the upper critical dimension 
calls upon a better understanding of their relation. 

On a practical level, understanding the FSC is an 
essential ingredient to infer correct equilibrium proper- 
ties in the thermodynamic limit from simulations of (in- 
evitably) small system size N££ Similarly, controlling the 
finite-size behavior is crucial to prove the existence of av- 
erages in random graph and combinatorial optimization 
problems.— 



II. FINITE SIZE CORRECTIONS IN SPIN 
GLASSES 

In this study we focus on the paradigmatic EA Ising 
spin glass defined by the Hamiltonian 2 . 

H = -^2 J l . ] cr l a J1 (1) 

(i,3) 

with spin variables Oj 6 {±1} coupled in a hyper-cubic 
lattice of size N = L d via nearest-neighbor bonds J,.j, 
randomly drawn from some distribution P ( J) of zero 
mean and variance (J 2 ) = Jo- The difficulty of comput- 
ing with any statistical accuracy the FSC of H at T = 
traces mainly to two disorder-specific complications: (1) 
Averages over many instances have to be taken, and (2) 
finding zero- or low-temperature states for each instance 
with certainty requires algorithms of exponential com- 
plexity. Thus, even with good heuristic algorithms^ at- 
tainable systems sizes are typically rather limited such 
that the form of corrections often must be guessedj 32 i 35 
Sometimes, theoretical arguments, such as those valid at 
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T = T c for SK^ can be used to justify the extrapolation 
of data. But there are no predictions yet for FSC below 
T c at mean-field level. For EA the conjectured scaling 
relation at T = (e.g., see Ref. [37j) is 



(2) 



For square lattices at T = = T c this was shown in 
Ref. [23| . Eq. ([2]) relies on the fundamental exponents**^ 
y(= 9) governing the "stiffness" of domain walls in low- 
temperature excitations. It is far more complicated to 
verify this relation in the glassy phase at T = for sys- 
tems with T c > 0, such as EA in d > 3M- The smallness 
of the stiffness exponent y in d = 3 obscures the distinc- 
tion between bulk- and domain-wall effects in the FSC. 
While y/d is found to increase in higher dimensions j2I 
computational limits on the number of degrees of free- 
dom, and hence on the length L, in higher-d simulations 
diminishes any such advantage. 

Here, we utilize diluted lattices to attain sufficiently 
large lengths L, in particular, in larger d to access the 
optimal scenario to distinguish scaling behaviors. Using 
the Extremal Optimization (EO) heuristi o 34 ! 39 to sam- 
ple many instances at T = 0, we find evidence in dimen- 
sions d = 4, . . . , 7 that the scaling behavior, controlled 
by domain-wall excitations, also applies when T < T C £Z- 



A. Cost versus Energy 

In ground-states simulations of finite-degree systems 
(like EA or BL) it is convenient to measure the "cost" 
Co, i.e., the absolute sum of unsatisfied bond- weights. 
The Hamiltonian for EA in Eq. ((T|) is a sum over all 
bonds, where this cost C contributes positively to the 
energy and satisfied bonds provide a sum S that con- 
tributes negatively, E = C — S. The absolute weight of 
all bonds, B = C + S is easily counted for each instance, 
amounting to a simple relation between cost and energy, 



E = -B + 2C 



(3) 



Averaging over undiluted (p — 1) lattices with discrete 
± Jo-bonds, the absolute weight of all bonds is fixed, B = 
dNJo. and E and C fluctuate identically. In contrast, 
for randomly diluted and/or a continuously distributed 
bonds, instances fluctuate normally around the ensemble 

mean, B ~ (B) + O (VF) with (B) = pdN (\J\) here. 
When only a small fraction of bonds contribute to C, 
as in these dilute systems, those trivial fluctuations in 
B contribute a statistical error to the measurement of 
E but not of C, although both exhibit the same FSC^i 
Thus, we prefer to extrapolate for the cost density 
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and use Eq. ^ with the energy density in the thermo- 
dynamic limit, (eo)^ = 2 (cq)^ — (B) /N, to obtain 
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When Co is intensive, in fact, only a measurement of 
C can provide the relevant FSC for E: Now (Co) N ~ 
(CoU + JjTTV 1 "" leads to Eq. © with (eo)^ - 
(— (B) + 2 (Co)oo) /N. A direct measurement of (eo) 
would merely yield the trivial constant — (B) /N with 
2(Cq) oc /N as equally trivial "correction". Note that 
any such finite-degree systems is very different from 
SK, for which (B) - JoN 2 /2 + 0(N) and (C ) ~ 

N 2 /4 (j - 2e SK /VN) so that by Eq. © e SK ~ 

(E) /iV 3 / 2 remains as ground-state energy density. 



B. Finite-Size Corrections from Domain Wall 

It is well-known that at low temperatures thermal 
excitations in Ising spin systems take the form of 
"droplets" which at size s require a fluctuation in the 
energy ~ s y , due to the build-up of interfacial energy. 
For y > 0, large droplets are inhibited, indicating the ex- 
istence of an ordered state, i.e., a phase transition must 
occur at some T c > 0. For instance, in an Ising ferromag- 
net, droplets must be compact to minimize interfacial en- 
ergies at s d ^ 1 , i.e., y = d — I, making di = 1 the lower 
critical dimension. In EA, such droplets are shaped irreg- 
ularly to take advantage of the bond-disorder, bounding 
y <(d- l)/2^ and in fact, it is d ; = 2.5^SL2S 

For EA at any bond-density p, the ground state 
energies are extensive on a finite-dimensional lattice, 
Eq (L) ~ L d = N. With periodic boundary conditions^ 
the FSC to this behavior is believed to arise primarily 
from frustration imposed by the bond disorder that mani- 
fests itself via droplet-like interfaces spanning the system. 
(An example are the "defect lines" connecting frustrated 
plaquettes in 2d.) Hence, 
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Table I: Stiffness exponents for Edwards- Anderson spin 
glasses for dimensions d — 3, ... ,7 obtained numerically from 
domain-wall excitations of ground states. The exponent yp 
refers to lattices at the bond-percolation threshold p c r~ an d y 
characterizes the T — fixed point in the glassy phase (since 
T c > for d > 3) for bond densities p c < p < 1.^2 The last 
column contains the measured values, denoted as Co, from the 
fit in Fig. H 
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-1.289(6) 
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0.920(4) 
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-1.84(2) 
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-2.01(4) 
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1.1(1) 
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-2.28(6) 
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0.823(7) 


0.91(5) 
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Figure 1: Plot of the approximate ground state cost per spin, 
(co), at bond-percolation (p = p c ) in dimension d — 3 (top) 
and 4 (bottom) on lattices up to L < 100. The data is plotted 
with respect to the presumed FSC, 1/L d " with u = 1 — yp/d, 
using y P = -1.29(1) in d = 3 and y P = -1.574(6) in d = 4, 
see Tab.|T] Data for L < 7 (way off scale to be shown) deviates 
from the asymptotic scaling. 



Rewriting Eq. © in terms of the energy density (eo) = 
(E ) /L d then yields Eq. j5J) with the FSC exponent w 
given by Eq. {2]), which is derived easily below dj with 
T c = 0, i.e., y < 0. 

Extending Eq. © naively above d; for < T < T c , 
i.e., y > 0, has not been tested with sufficient data even 
in d = 3 i 37 i 40 ' Since finding ground states becomes a 
hard combinatorial problem requiring heuristics, attain- 
able system sizes L have been too small to savely distin- 
guish domain-wall generated from higher-order correc- 
tions such as bulk effects (u = 1), since y is small. E.g., 
in d — 3, the lowest dimension above di, for which the 
largest systems sizes (L < 12) are accessible, we have 
y = 0.24(1), i.e., y/d « 0.08 and cj = 1 - y/d « 0.92. 

Using bond-diluted lattices, we have previously deter- 
mined accurate values for y from the domain- wall energy 
in d = 3, . . . , 7i 19 ' 20 For one, dilute lattices provide large 
length scales L and higher d for which ground states can 
still be obtained in good approximation. For our pur- 
poses, a higher d also offers a more distinct ratio for 
y/d, see Tab. fl] Already in d = 4 it is y = 0.61(1), i.e., 
y/A as 0.45, about double the value in d = 3. 



III. EDWARDS-ANDERSON MODEL AT THE 
PERCOLATION THRESHOLD 

As a demonstration of our approach, we first treat the 
problem on strongly diluted EA-lattices exactly at p = 
p c , previously studied in Ref. [22|]. The fractal lattice at 
p c is too ramified to sustain order and y = yp < for 
all d42 Polynomial-time algorithms have been used to 
achieve system sizes up to N « 10 8 , i.e., L < 11 in d = 7, 
and the results for yp are recounted in Tab.fl] With the 
available data, we focus on demonstrating its consistency 
with Eq. ^ by plotting it on the appropriate scale. The 
results for the extrapolation of the ground state cost for 
increasing system size L according to Eq. Q are shown 
in Fig. []] for d = 3 and 4. Indeed, the linearity of the 
plot on the scale 1/L du) with u = 1 — yp/d, as given 
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Figure 2: Plot of the average ground state cost (Co) as a func- 
tion of 1/N 1 / 3 , obtained with the reduction method described 
in Refs. [l9ll2(j |. for percolating random graphs Gjv iP (degree 
Np = 1) of sizes up to iV = 2 19 , using a continuous (Gaus- 
sian) bond distribution. Apparently, the total cost itself is 
intensive and approached with (Co) N ~ (Co)^ + |T /N 1 ^. 
Inset: Corresponding plot for discrete bonds, exhibiting en- 
tirely distinct scaling, (C ) N ~ (Co}^ ln(JV) + T/2. 



in Tab. U for p = p c , supports Eq. $Z§, as expected for 
T c = and yp < 0^- However, it is interesting to note 
that the data in Tab. U suggests w — > 4/3 for d — s- oo, 
as depicted in Fig. [3] In the large-d limit, a lattice at 
p = p c ~ 1 / (2d) is similar to an ordinary random graph^ 
Gjv, p at percolation, Np — > 1^ this would predict FSC 
of the form 1/iV 4 / 3 for the cost density of percolating 
Gjv.p- As Fig. [2]shows, we do find an FSC with uj = 4/3 
when using Gaussian bonds, except for the fact that the 
cost (Co) itself in such a limit has become intensive. The 
energy remains trivially extensive, Eq ^ — B 3> Co, as 
discussed above. 

Fig. [5] also shows the corresponding result using a dis- 
crete (±Jo) bond distribution but the same graph en- 
semble. Note that the change in distribution affects an 
entirely different scaling. The leading-order difference 
can be explained by envisioning Gn,p near percolation as 
a tree possessing O (In N) weakly interlinked loops, each 
of length O (In A), with half of those (ld-like) loops frus- 
trated in a single bond: The total cost for discrete bonds 
oi fixed weight |J| = Jo is (Co) N ~ Jo In A, while for 
continuous bonds this cost is mitigated by selecting the 
weakest weight, |J| ~ Jo/ In A, in the loop and (Co) N 
approaches a constant value. 

We conclude for p = p c that our results in finite d 
connect well with those for d — > oo, as long as we choose 
a continuous bond distribution in the mean-field limit. 
While non-universality is not unexpected when T c = O^ 3 - 
we will encounter it also for FSC when T c > because 
FSC do not affect the thermodynamic limit. 
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Figure 3: Plot of 1 — y/d and the measured values, Q, in Tab.|l] 
as a function of inverse dimension 1/d. The vertical lines 
marks the upper and lower critical dimensions of EA at d u — 6 
(left) and di ~ 5/2 (right), at which y — 0. Possible mean- 
field predictions to hold above d u are marked by horizontal 
lines: u) = 4/3 for lattices at p = p c , u = 5/6 for EA- lattices 
at p > p c , and w = 2/3 as predicted from simulations of SK. 
The extrapolations in Figs. Q] and U suggest the validity of 
Eq. ((2)1 for all d, and it seems difficult to relate the FSC in 
SK to those of EA. 



IV. EDWARDS-ANDERSON MODEL IN THE 
GLASSY STATE 

Parallel to the above case, we analyze our data for di- 
lute EA-lattices in the glassy regime, p > p c . However, 
this data had to be obtained in far more laborious simu- 
lations. While starting from lattices of comparable size, 
up to N ~ 10 8 , the methods from Ref. [22| merely reduce 
each instance to a remainder of about Nr < 10 3 variables 
at the chosen values of p > p c . (The choice of p is crucial: 
larger values dramatically increase the remainder sizes, 
smaller values create long transients due to the proximity 
of p c .) The cost Co of each remainder had to be optimized 
by local search with EO^ run with O (N R ) updates^ 
Depending on system size L and dimension d, disorder 
averaging required 10 4 — 10 7 instances to reduce statis- 
tical errors. In contrast to our earlier simulations ; 1 ^ 20 
far more statistical accuracy for each energy average and 
new implementations to handle very large lattices effi- 
ciently are required here4^ While there is never a guar- 
antee for finding actual ground states, we have applied 
similar safeguards to hold systematic errors below sta- 
tistical errors as in previous investigations! 26 ' 2 ^ 34 ' 39 For 
instance, in d = 7 at L = 11 we have tested subsets of 
100 instances at run-times of twice the normal setting 
and did not find a single deviation. 

To avoid the above-mentioned problem of wide separa- 
tion between data points at L and L + 1 in higher d, we 
have also undertaken to simulate asymmetric instances 
of size N = L d - m (L + l) m for m = 1, 2, . . . , d- 1. While 
we display that data as well, it proved unsuitable to fit. 



It merely serves to demonstrate that even such a minor 
variation in rather large, random instances resulted in 
sufficiently large systematic distortions from the m = 
data to be detectable with our heuristic. 

The results of this effort for d = 3, . . . , 7 are presented 
in Fig. 2J In each case, the raw data in the asymptotic 
scaling regime has been fit to the form in Eq. (fj| to yield 
w = 0.914(5) and (co)^ « 0.0425063 with x 2 /ndf « 1.3 
(d = 3), Q = 0.82(1) and (co)^ « 0.0063720 with 
X 2 /ndf « 0.77 (d = 4), Q = 0.81(1) and (cq)^ « 
0.0022543 with x 2 /ndf « 1.8 (d = 5), u> = 0.82(2) 
and (co)^ « 0.0014851 with x 2 /ndf w 0.8 (d = 6), 
u) = 0.91(5) and (co)^ w 0.00061055 (d = 7). Each 
fitted value of w is close to 1 — y/d, see Tab. U and Fig. [3l 
except for d = 7. For the largest d, the value of tu trends 
higher, although it is difficult to judge whether asymp- 
totic behavior has been reached. Fixing u> = 1 — y/d 
and using a l/L da; -scale as in Fig. [TJ a two-parameter 
linear fit (not shown) yields virtually identical values: 
(co)^ « 0.0063735 (d = 4), (co)^ m 0.0022546 (d = 5), 
(co)^ « 0.0014852 (d = 6), and (co)^ « 0.00061047 
(d = 7) . In Fig. 2] we have subtracted the fitetd value of 
( c o)oo f rom each data point to be able to plot the remain- 
ing scaling form on a double-logarithmic plot, for better 
visibility. 

To assess the effect of higher order corrections on our 
fits, we fixed the parameters of the l st -order fitSS and 
added another scaling correction term to Eq. ^ for a 
2 nd -order fit, now including all points down to the small- 
est size, L = 3. In each case, these corrections scaled 
with an exponent of w 1.3 — 1.6, about twice that of the 
leading correction. Shown as dotted lines in Fig. 0] these 
higher-order corrections have virtually no impact on the 
asymptotic behavior. We have tested that the existence 
of such a rapidly vanishing correction does not affect our 
results by more than 2% in Co. 



V. CONCLUSIONS 

We conclude, first and foremost, that the ground-state 
FSC in EA are consistent with Eq. for all d < d u = 6, 
whether T c = (p = p c ) or T c > (p > p c ). The con- 
nection with mean-field predictions is more complicated. 
The FSC exponent w for EA does not approach its equiv- 
alent of lusk » 2/3 in SK 2 ^ 2 ^ at d u , by far. As Fig. [3] 
shows, its mean-field limit seems more aligned with the 
result of u> w 0.8 found for sparse-degree BL having a con- 
tinuous bond distribution^ There is no obvious reason 
why those BL would provide a more appropriate mean- 
field limit for EA, except for their sparse degree. In the 
computationally simpler case of EA at p = p c , a simi- 
lar approach for FSC to a sparse mean-field limit is even 
more apparent. 

Beyond FSC, ground-state energy fluctuations provide 
another example where its exponent p in EA does not 
reach the corresponding SK value for large d. It has 
been proven 4 ^ that those fluctuations are normal (i.e., 
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Figure 4: Double-logarithmic plots of the approximate ground state cost per spin, (co) L = (Co) L /L d , reduced by its numerically 
extrapolated limit value (co}^ = lmiL-^oo {co) L , versus l/L d for dimensions d = 3, . . . , 7. Subtracting off (co)^ highlights the 
FSC-term in Eq. ©. Simulated lattice sizes range from L = 3 to L < 25 (d = 3), L < 20 (d = 4), L < 15 (d = 5), and L < 10 
(d = 6, 7) at respective bond densities p = 0.185 (d = 4), p = 0.13 (d = 5), p = 0.105 (d = 6), and p = 0.082 (d = 7). The 
l st -order fit (solid line) to obtain (co) and to using Eq. Q only utilized the fitted data-points. The measured values for Q 
are reported in Tab. [I] The 2 nd -order fit (dotted line), involving all data for d = 4, . . . , 7, held previously fitted parameters 
fixed but added a further correction term to Eq. Q. This correction shows little effect on the asymptotic scaling, as solid and 
dotted lines virtually coincide. The dashed line for any "bulk" corrections of order 1/L d is provided for reference only. Data 
( ) in d = 5, 6 obtained from asymmetric instances (labeled "asym") was exclude from the fits; it is displayed to demonstrate 
its consistency, being systematically "arced" above the symmetric data. 



p = 1/2) for all d in EA, whereas in SK it has a value of 
p = 5/6. 45 (In fact, it has been argued 4 ^ that p = 1 — y/d 
for d > d u , for which Eq. ((2|) implies co = p = 5/6 at d u , 
consistent with our data in Tab.|I]and Fig. [3]) Again, BL 
with Gaussian bonds appear to have the right mean-field 
limit for EA, exhibiting p = 1/2M- 

To sort out these behaviors we first note that both, 
lo and p, merely describe properties of the asymptotic 
approach to the thermodynamic limit and, hence, may 
exhibit non- universal (i.e., distribution-dependent) be- 
havior, even though they can be related to other, uni- 
versal exponents such as y. Non-universal, distribution- 
dependent behavior is far more prevalent for sparse 
mean-field spin glasses than for SKj 27 i 31 For such prop- 



erties, it ought not surprise that these sparse systems 
provide a better high-dimensional limit for EA than 
SK. After all, in studying EA we first take the limit 
N = L d — > oo at fixed d (for a sparse degree of 2d) be- 
fore we consider d — > oo, whereas SK corresponds to the 
correlated limit, N — > oo, d — > oo with degree 2d ~ N 
held fixed. In disordered systems, the order in which 
those limits are taken might well matter. However, the 
quest to understand the extend or possible break-down of 
RSB, which is known to describe the glassy state in sparse 
mean-field model a 24 ' 25 as well as SK, remains undimin- 
ished. 

This work has been supported by grant DMR-0812204 
from the National Science Foundation. 
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